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A nondegenerate two-photon Jaynes-Cummings model is investigated where the leakage of photon 
through the cavity is taken into account. The effect of cavity damping on the mean photon number, 
atomic populations, field statistics and both field and atomic squeezing is considered on the basis 
of master equation in dressed-state approximation for initial coherent fields and excited atom. 
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I. INTRODUCTION 

ly.^ ■ Over the last two decades much attention has been focused on the properties of the dissipative variants of the 
^vq Jaynes-Cummings model (JCM). The theoretical efforts has been stimulated by experimental progress in investigation 
of the interaction of a single atom with electromagnetic field inside the cavity [lj . The experiments with highly excited 
t— I ' Rydberg atoms allowed some of the predictions of the extended version of JCM to be proved. Besides the experimental 
drive, there also exists a theoretical motivation to include relevant damping mechanism to JCM because its dynamics 
becomes more interesting. The dissipative effects in JCM caused by the e nerg y exchange between the system and 
environment have been studied both analytically 0-0 and numerically 0-0I- Last few years the JCM with phase 
damping, as applied to decoherence and entanglement, has been also treated intensively [ill H^. 

It's known that two-photon processes are very important in atomic systems due to high degree of correlation 
. between the emitted photons. Hence, one valuable extension of the JCM is well-known two-photon JCM. With the 
experimental realization of two-photon micromaser 0] the dissipative two-photon JCM has attracted a great deal 
of attention ^ is worth nothing that in all mentioned works the dissipative dynamics of the degenerate 

£v' two-photon JCM has been under consideration. In order to advance one step further in the investigation of two- 
photon processes, the JCM with two-mode two-photon interaction or nondegenerate two-photon JCM (NTPJCM) 
were proposed. A remarkable feature of such a model is that one mode can be used to affect on the other mode. 
The lossless NTPJCM has been used to study time evolution of the atomic and photon operators, the second-order 
coherence function, the one- and two-mode squeezing, the atomic dipole squeezing, the emission spectra and quantum 
entropy and entanglement without and with consideration of Stark-shifts 0]-[3l]]. The influence of phase damping on 
£> . nonclassical properties of NTPJCM has been considered in • The effect of a cavity damping on the time behaviour 
of the atomic population in the special case when the fields are initially in the two-mode squeezed vacuum has been 
taken investigated by Gou [l9|. It's of great interest to investigate the role of energy dissipation in dynamics of the 
NTPJCM for arbitrary fields states. 

II. MODEL HAMILTONIAN AND KINETIC EQUATIONS 

The nondegenerate two-photon Jaynes-Cummings model is an effective two-level atom with upper and lover states 
denoted |e) and \g) respectively interacting with two modes of quantum electromagnetic field with frequencies uj\ and 
u>2 through two-photon transition. The Hamiltonian for such a system in dipole and RWA approximation is 

H = TiloqR z + hiuj-iOi a\ + ui20'2 a ^) + hg{a\a2R + + a^a^ R~), (1) 

where u>o is the atomic transition frequency, uii and 0J2 are the cavity mode frequencies, af (ai) are the creation 
(annihilation) operator of the photon (i = 1,2), R z is the inversion population operator, are the operators 
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describing the transitions between the upper and lower levels and g is the atom- field coupling constant. We have 
ignored the Stark shift caused by the intermediate level and denoted the detuning parameter as 

which satisfy the condition A -C u>o, k>2- 

In order to describe dissipation one has to treat the system as open. In this paper we take into account only the 
field mode damping and ignore the atom damping. The cavity is assumed to be at zero temperature. Then, the 
master equation for the density matrix of combined (atom-field) system is 



^ = -ij h [H, p] - kj (afaip - 2 mp af + pajai) , (2) 

i=l 

where 2fcj (i = 1, 2) are the rates of photon leakage from the cavity. For the sake of simplicity we put ki = fc 2 = k. 
Using the representation 

W(t) = e* m p{t)e-k H \ 0(t) = e^ Ht Oe~^ m , 
where O is an arbitrary operator of combined system, one can rewrite the master equation (2) in the form 

^ = - Y, k ( a t^W - 2 aiW af + Wa+cn) , (3) 

i— 1 

To solve Eq.(3) we have used the so-called dressed-states representation, i.e. representation consisting of the 
complete set of hamiltonian eigenstates. For lossless cavity the full set of dressed states are 

= ^|+, ni ,n 2 ) ± ^ |-,m + l,n 2 + 1), (4) 

with eigenvalues 

E n fo4*n\ri2 f- ^nin2 j 

where 

(Amna = + \) + ^2(^2 + i), 



/A 2 

V nin2 = d — + ,g 2 K + l)(n 2 + 1), 8{n) = A/n(n). 

Here \a;n > refers to a state with n photons in the cavity field mode and the atom in the excited (a = +) or in the 
ground (a = — ) state 

\a; n >= \a >a \n >f, 

where n = 0, 1, 2, 

For finite-Q cavity the above states (4) should be added with the states: \^/\} = \— , 1,0), E — Huji — ^Hcuq, 
\¥ 2 ) = |-,0,1), E = Hlj 2 - ^Hlj , 

l^l) = |— , 0, 0), E = — ifiwo, which take into account the photon leakage with no atom change. 
Using the secular approximation which holds for 2fcn 2 <C g\Jn + 1 i.e. neglecting the oscillatory terms, the 
equations for the diagonal elements of density matrix W are found to be 

«„JW1<„ 2 ) = -* {2(m +n 2 ) +7^ l2 «„J^K„ 2 )- (5) 



9 2 {n 2 + 1) 



Til + 1 7ni + l,n 2 



ni 



Tnin 2 



^ni+l,na 7 ni + l,n 2 



\ ni+l,na 
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9 2 (n 2 + 1) 



m + 1 7„ 1+ i,n 2 "1 + 2 7if lTl2 



^m+l,n 2 7,Ti + l,n 2 



/\[r + 



W^m + l.n,}- 



3 2 (ni + l) 



n 2 + l7 ni ,„ 2+ i n 2 + 2 7^ in2 



^niH2 7nii2 



^ni,n 2 +l 7,ti,n 2 +l 



(*n 1> n 2+ l|W|*n 1 ,n 2+ l)- 



<? 2 (m + i) 



n 2 + l7n 1 ,n 2 + l _ "2 + 2 7^ 



7"l«2 



^m,n 2 +l 7ni ,n 2 +l 



(*^,„ a + l|W|*^,n a + l>. 



where 



7mn 2 



1± 



2fL 



and 



"oo7oo 



(6) 



7^2<*ool^l*oo) 
"oo7oo 



(¥ 3 \w\¥ 3 ) = 2k((¥ 1 \w\¥ 1 



(i = l,2), 

* 2 |^l*2»- 



(7) 



The equations for off-diagonal elements of IF with nonzero right-hand sides are 

(*$ ina \W{t)\9* na ) = -2k(n 1+ n 2 + l)^t in2 \W\^ in2 ). 

The solutions of equations (8) are 

(^ in JW(t)|^ in2 ) = (<„JW(0)|^ in2 )exp{-2fci(n 1+ n 2 + l)} 

and the solutions of equations (5)- (7) may be obtained only numerically. For this purpose one can assumed that initially 
there is an upper limit on the number of photons N\ and N 2 in both of cavity modes so that (^nm 2 \^(^)\^nm 2 ) = 
for m > N\,n 2 > N 2 . This implies that these matrix elements are zero for all t since the cavity cannot add to the 
photon numbers. Then, one can start with n\ — N\ + 1 and n 2 = N 2 + 1 and iterate equations (5), (6) for smaller 
values of photon numbers untill n\ = n 2 = 0. If there is no upper limit on the initial numbers of photons in the 
system the numbers N\ and N 2 must be taken large enough for the mean values of observables to calculate with the 
appropriate accuracy. These quantities may be obtained in the standard manner 

(0(t)) =S P 0{t)W{t). 

The solutions of the Eqs. (5)-(7) for arbitrary initial states of atom and field can result from numerical calculations. 
We consider below the NTPJCM with the atom initially in the excite state and the fields in coherent states. 



III. RESULTS AND DISCUSSIONS FOR COHERENT INPUT 

The initial density matrix W(0) for atom in the excited state and the fields in the coherent states is 
W(0) = P^Emf^ 2 m+ \m+ i + 7 - 2 \m- i+ 

n 1 n 2 \ + 7n 1 n 2 7n 1 n 2 l^'n 1 n 2 ) (^'ti 1 ti 2 I) ' 

where 



n 

Pn i =exp(-n i )- L - (i = 1,2). 
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First consider the time behaviour of mean photon numbers and mean atomic populations 

OO OO ^ 

W(*)> = E E { -7i- 1 „ a 7^ ina «„ a |W|*- na )c 0S (2n nina t)+ 



+ (n, + (*^ ina |W|*- n2 ) + (n, + (*+ „ 2 |W|*+„ 2 ) J + 

+ = 1,2), 

OO OO ( _|_ 2 — 2 

m— 0ri2— L 



+ 7i„ 2 7n 1 „ a (*n 1 „ a l^l*n 1 „ a )cOs(2fi„ 1 „ 2 i) J 



OO OO 



<*»(*)> = E E { ^(<nJW|< n J + ^p(*- n2 |^|*- n2 )- 



ni- ri2— 



z=l 



The mean populations of the excited atomic state in the presence of the two modes of coherent state are plotted 
in Figs. 1 - 4 for various values of (N\), (N 2 ), S and k. For small values of 5 and k the phenomena of quantum 
revivals and quantum collapses of the Rabi oscillation appear. They are not as regular as those in onc-photon or 
degenerate two-photon case. It can be seen that the amplitudes of the revival oscillations decrease as a consequence 
of the cavity damping and detuning. In the case of strong damping, the cavity losses are so large that no collapses 
or revivals phenomena may appear. As a result of the computer simulations it can be said that the decay time and 
serene duration time for atomic populations are directly affected by the initial photon numbers in the cavity modes. 
The serene duration time decreases as (N\), (N2) decreases and mean populations manifests the more fluctuating 
behaviour. For large field intensities the detuning influences both revival amplitudes and serene duration time, as 
well as quasi-stationary atomic population value. 

Since two photons are absorbed and/or emitted by the atom simultaneously in the cavity, one can tell that the 
behavior of photon numbers of mode 2 shows the exactly same manner as for the mode 1. Therefore in Figs. 5 -8 we 
have plotted the mean photon number for first cavity mode. The mean photon numbers exhibit the same pattern 
of collapse and revival as the atomic population. A comparison of Figs. 1-8 shows that the photon numbers more 
significantly affected by cavity damping than the atomic population. In the case of strong damping or detuning the 
mean photon number decays exponentially. 

Perhaps a better appreciation of the statistics can be had by examining second-order correlation function as 
a function of time t. The second-order correlation functions for two cavity fields may be defined as 

G (2) {t)= ((^(t)) 2 a 2 (t))-(a+(t)a(tr 



(a+(t)a(t))* 

For strictly coherent field G^ 2 \0) = whereas negative values of lead to the antibunching of the field. In 
dressed-state representation the second-order correlation function becomes 

1 ~ ~ (( o A 



+ (^+71!— — W in2 |W|*+ ina >- 

- 2n i 7+ na 7- n2 <*+ n2 |W|*- n2 )c£w(2n nin2 t) \ - 1. 



(2) 

In Figs. 9 and 10 we plot G\ for different detuning and damping parameters and large initial field intensities. For 

undamped resonant cavity collapses and revivals appear. In the case of collapse in the absence of detuning (or cavity 
(2) 

damping), G\ < , the oscillations show both bunching and antibunching features. In the cases of nonzero damping 
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and detuning the antibunching effects almost disappear. The amplitudes of the revivals oscillations will be damped 
by small cavity losses to the extent that antibunching appears only at the very beginning of the time evolution. For 

(2) 

large detuning the second-order correlation function G\ > for every t. 

Finally, we study the field and atomic squeezing. In order to investigate the squeezing properties of the radiation 
field we define the slowly varying Hcrmitian quadrature operators for fields 

= ^(a i e^ t + afe-^ t ), 
(i = l,2). 

The commutation of x[ l} and X^ is [x[ l \x^] = i/2. The variances {Axf) 2 = {{xf ] ) 2 ) - (xf) 2 (j = 1,2) 
satisfy the uncertainty relation (AX^^^AX^) 2 > 1/16. For the vacuum and coherent states of the field the 



variances are equal 1/4. The field is in a squeezed state if there takes place (AX- ! ') 2 < 1/4 for cither j = 1 or 2. 



The condition for squeezing in the jth quadrature AX can be written simply as 



Sf < 1, 



where the squeezing factor is 



Sf=4AXj <) (i = 1,2). 



For the sake of definiteness we study the squeezing properties of the first cavity mode. In terms of the photon 
operators, one can readily find that the squeezing parameter of the first quadrature component and for the first cavity 
mode may be written as 

S = (a 2 ) + (af) + 2(af ai ) - «a+ + fll )) 2 + 1, 

where 

oo 

(ai) + (af) = + 2 {(j+ in2 Jn 1+ 2n 2 Vm + 1 + 7^ in2 7ni+2n 2 >/ni + 3)x 

ni ,n2— 

X (^n 1 +2n 2 \ W \^n 1 n 2 ) cos (Pn 1 +2n 2 

+ (T^ n 2 7n, +2 n 2 + 1 - 7n, n 2 7ni +2 n 2 V^l + 3) X 

X {*n 1+ 2n 2 \W\*Z in2 }cOs{[n ni+2 „ 2 + 0„ in2 ]i) + 

+ (7nj n 2 7nj +2 n 2 V™1 + 1 + In, n 2 In, +2 n 2 V^l + 3) X 

+ (7 m „ 2 7nj +2 n 2 V^l + 1 - 7+ „ 2 7~ +2 „ 2 V^l + 3) X 
X (*i 1+ 2n 2 \W\*- in2 )cOs{[n ni+2n2 + fi ni „ 2 ]t) }, 

OO 

( ai ) + (a+) = ^ {(7n 1 n 2 7n 1 + ln 2 + 1 + 7„ 1 n 2 7,7 1 + l n 2 + 2)x 

"1 ,77.2=0 

+ (7n 1 n 2 7n 1 + l n 2 V"l + 1 - 7nin 2 7n 1 +l n 2 V"l + 2) X 

x Jw|*+„ 2 )c 0S ([Q„ 1+1 „ 2 + n nin2 ]t) + 

+ (7n 1 n 2 7rT 1 +ln 2 Vni + 1 + 7rt 1 n 2 7,t 1 + l „ 2 V^l + 2) X 

x Jw|*-„ 2 )cos([Q„ 1+1 „ 2 - fi nin2 ]t) + 

+ (7n 1 n 2 7^" 1 + l„ 2 V»l + 1 - 7n in2 7ni + ln 2 V™1 + 3) X 

x (*^ 1+ i„ 2 |w|*- in2 )cos([n ni+ i„ 2 +Q nin2 ]t) }, 
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and (a+ai) = {Nx(t)). 

Values S(t) < 1 imply squeezing in the first quadrature component. In Figs. 11 and 12 we have plotted S(t) for 
different values of k and 5 and large input field intensities. For k — there is some amount of squeezing which appears 
only at the very beginning of the time and this decreases with increasing of 6 in contrast with degenerate two-photon 
JCM. But this model does not give rise to significant squeezing. For nonidcal cavity the effect is vanished. So, the 
cavity damping is seen to have an appreciable effect on the squeezing properties. 

To investigate atomic squeezing we introduce the dispersive and absorptive components of slowly varying atomic 
dipole moment |32( 



and 



o- 2 = JL(J2+ e -«"»* - re" 9 '), 



respectively. They obey the commutation relation 



[01,02] = -03 



and the corresponding uncertainty relation 



(A ( r 1 ) 2 (A<7 2 ) 2 > i(o 3 ) 2 . 



The atomic state is said to be squeezed when <j\ or er 2 satisfies the relation 



4 

Since 



(Ao-0 3 < 7IWI, (i = l,2) (9) 



(A CTl ) 2 = i - (Jfe<a)e«**) 2 , 



(A^) 2 = I - {Im{a)e^)\ 



the condition described by Eq. (9) can be rewritten as 

1 - A(Re(a)e-^ "*) 2 

Fl= <X 

or 

^2 - H Tj < 1- 

IWI 

for squeezing in the dispersive or absorptive component of the dipole moment. 
Here 

iAt 00 

( a ) e -^0t = £ \^+ ,~- e -**(°"l + l "2 + 1-^1 "2 W+ |^|*+ n ^) 

\ / 2 / ■ > (_ 'ni + lna + l lni"2 \^niri 2 l" I m + 1 n 2 + l/ 

n± ,ri2— 

_ 7 + 7 + -^(n„ 1+1 „ 2+1 +n, n „ 2 )/^- 1^1^+ \ + 

'ni + l n2 + l 'ni 712 \ ni n2 1 m+lna-l-l/ 1 

— j— 7" ,. ,,-y- e «*(f2"i + i"2 + i+ fi "i"2)($+ n 
' Ini+lri2 + l mi n 2 \ I ni n 2 l " I ^ ni+1 n 2 + l' 

- 7~, n ^ n 7 + e It ( fi "i+i"2+i-""i"2)(^- IVFI^ 1 

/ni+lri2 + l 'riin 2 \ x iVinjl" I 711 + I712+I/ | ' 

and (173) = (i? e ) — (Rg) is atomic inversion in terms of dressed states representation. 
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The results for the squeezing parameters for different values of k and 6 and moderate input intensities ((rii) — 
15, (712) = 10) have been shown on Figs. 13-16. The dispersive component F% does not squeeze at the very beginning 
of the time, the absorptive component F2, on the other hand, goes below 1 with virtually no time delay. 

Both F\ and Fi shows squeezing recurrently only for small times of the atom-field interaction. As that has been 
mentioned earlier |23j that squeezing does not show up in either case until after (iVj) > 7.0. The amount of squeezing 
decreases with increase of damping parameter k. If the parameter k is large enough the squeezing in dispersive 
component vanishes and the squeezing in absorptive component of dipole moment undergoing only one squeezing 
minimum. For high input intensities (not shown in Graphs) the influence of damping parameter on squeezing amount 
is more dramatic than one for medium input intensities. And with increasing the detuning parameter S the amount 
of squeezing in dispersive and absorptive component also decreases while the time interval, for which the squeezing 
appears, increases. 

IV. SUMMARY 

In this paper we have investigated the nondegenerate two-photon Jaynes-Cummings model with damping and 
detuning. We have set kinetic equations for density matrix of the considered system in secular approximation and 
with using the dressed-state representation. These equations are solved numerically for different detuning of atomic 
levels and damping parameter. On the basis of these equation the analysis of the dynamical behaviour of mean 
values of atomic populations, mean photon numbers, field coherence and squeezing has been carried out. The effects 
of cavity damping will significantly attenuate the amplitudes of mean atomic populations, mean photon number 
revivals. The revivals of the second-order photon correlation function will be damped by small cavity losses to the 
extent than antibunching appears only at the very beginning of the time. Cavity damping has an appreciable effect 
on the squeezing properties of fields and atomic dipole moment. For moderate damping parameter the field squeezing 
is vanished and the amount of atomic squeezing sharply decreases. In this paper we shall restrict our consideration 
to the coherent input for fields and ignore the Stark shift. A further discussion on dissipative NTPJCM including the 
consideration of initial squeezed and thermal states for cavity fields and Stark shift is planned to be reported in the 
subsequent paper. 
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FIG. 2: Atomic population of an excited level for (JVi) = (N2) = 5, 8 — 10 and 1. fc = 0, 2. fc = 
0.001, 3. k = 0.01. Curve 1 corresponds to value (R e (t)) + 0.2, and curve 3 corresponds to {R?(t)) — 0.2 
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FIG. 3: Atomic population of an excited level for {N^ = (N 2 ) = 30, k = 0.001 and 1.5 = 0; 2. 5 = 20; 3. 5 = 100. 
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FIG. 4: Atomic population of an excited level for (Ni) = (N 2 ) = 30, 8 = 10 and 1. fe = 0; 2. fe = 0.0001; 3. k = 0.001; 4. k = 
0.01. Curve 1 corresponds to value (R e (t)) — 0.2, curve 3 corresponds to {R e (t)) + 0.2, curve 4 corresponds to { R e (t)) + 0.4. 
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FIG. 5: Mean photon number in the first field mode for {Ni} = (N 2 ) = 5, k = 0.001 and 1. 5 = 0; 2. 5 = 10; 3. 8 = 100. 
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FIG. 6: Mean photon number in the first field mode for (Ni) = (iV 2 > = 5, <5 = 10 and 1. k = 0, 2. fc = 0.001, 3. k = 0.01. 



30 „ 
30 

A 

3 30 
V 30 

































































1 


















8 














2 

\— 3— 



















2.5 5 7.5 10 12 .5 15 17.5 20 



FIG. 7: Mean photon number in the first field mode for (Ni) = (JV 2 ) = 30, k = 0.001 and 1. 6 = 0; 2. <5 = 20; 3. 5 = 100. 




10 



1:5 
t 



20 25 20 



FIG. 8: Mean photon number in the first field mode for (A^i) = {N2} = 
10 and 1. k = 0; 2. k = 0.0001; 3. k = 0.001; 4. fe 
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FIG. 11: Squeezing in the first field mode for (Ni) = (N 3 ) = 50, k = and 1.5= 10; 2. 5 = 100. 
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FIG. 12: Squeezing in the first field mode for (N x ) = {N 2 } = 50, 5 = 10 and 1. fc = 0; 2. fc = 0.0001. 
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FIG. 13: Atomic dipole moment dispersive component for (Ni) = 15, (N 2 ) = 10, k = 0.001 and 1. 5 = 50 ; 2.5 = 100. 




FIG. 14: Atomic dipole moment absorptive component for (Nx) = 15, (N 2 ) = 10, k = 0.001 and 1. 5 = 50; 2. 6 = 100. 



